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INTRODUCTION 


This  memorandum  is  a  review  and  summary  of  the  prominent  computational 
and  theoretical  approaches  to  the  modeling  of  reverberation  processes  in  the 
sea.  This  review  is  neither  exhaustive  nor  does  it  attend  with  the  same 
attention  to  every  topic  of  reverberation  attributes  and  characteristics.  It 
serves  as  a  base  line  of  documenting  the  work  in  reverberation  since  Fortuin's 
review  article1  and  the  more  didactic  works  preceding  Fortuin's  article. 
Another  aim  of  the  review  is  to  identify  the  essential  assumptions  of  the 
methodologies  currently  receiving  attention,  an  example  of  which  is  Wilson's 
work2  at  ARL  on  validating  Middleton's3  basic  theories  on  covariance  and 
statistical  properties  of  backscattering.  The  review  identifies  some  areas 
that  have  been  underutilized,  e.g.,  zero  and  level  crossing  characteristics, 
and  other  areas  that  have  benefitted  from  technological  development  (echo-echo 
correlation,  robust  detection).  Future  work  is  planned  to  review  in  more 
depth  promising  areas  and  to  identify  and  test  novel  concepts  of  character¬ 
izing  stochastic  fields  in  general  and  reverberation  in  particular. 

The  existence  of  reverberation,  in  underwater  sonar  applications,  is  one 
of  the  fundamental  limiting  factors  in  recognition  (detection  and  classifica¬ 
tion)  of  objects  one  may  wish  to  be  aware  of.  One  may  desire  information  as 
to  range,  bearing,  and  classification  of  objects  or  bottom  relief.  Similar 
infonaation  is  of  interest  in  radar  applications  where  the  infonaation  is 
carried  typically  by  electromagnetic  waves. 

Reverberation  arises  from  the  interaction  of  an  acoustic  or  an  electro¬ 
magnetic  energy  field,  ax>deled  by  waves,  with  boundaries  and  volume  scatterers. 

The  scattering  regime  is  typically  divided  into  the  geometric  and  dif¬ 
fractive  doaiains.  The  geometric  domain  is  characterized  by  features  much 
larger  than  the  wavelength  of  the  incident  radiation  while  the  diffractive 
doaiain's  features  are  of  the  same  order  of  magnitude  of  the  wavelength. 

Analysis  in  the  geometric  domain  is  typically  simplified  by  the  use  of 
ray  concepts.  This  is  of  great  use  in  describing  propagation  over  large 
distances . 

With  sound  waves  interacting  with  the  sea  surface,  the  incident  wave 
length  and  roughness  scale  are  often  of  the  same  order  of  magnitude  and  the 
interaction  is  characterized  by  diffraction.  This  may  also  be  true  of  bottom 
surface  interactions. 

There  are  two  boundary  conditions  that  idealize  the  extremes.  One  is  the 
free,  elastic,  or  pressure  release  boundary.  The  free  boundary  is  represented 
by  the  sea  surface.  The  second  is  the  rigid  boundary;  e.g.,  rock  bottom. 
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Thai  wave  velocity  potential  vanishes  on  a  free  boundary  while  the  direc¬ 
tional  derivative  of  the  potential  vanishes  on  the  rigid  boundary. 

The  interaction  of  sound  with  boundaries  depends  on: 

1.  Frequency  of  incident  wave 

2.  Geometry  of  source  and  receiver 

3.  Transient  nature  of  incident  wave 

4.  Temporal  variation  of  the  boundary.  The  sea  surface  varies 
in  time,  raindrops,  leaves  change  position  with  wind. 

5.  Geometry  of  boundary 

6.  Material  characteistics  of  boundaries. 

Reverberation  is  often  referred  to  as  unwanted  returned  energy  which 
■asks  the  wanted  signal  or  echo.  Urick4  calls  the  sua  of  all  the  scattering 
contributions  froa  all  the  scatterers  reverberation. 

Reverberation  has  several  features  which  distinguish  it  froa  noise.5 
Paramount  aaong  these  is  the  fact  that  reverberation  is  generated  by  the  sonar 
itself.  Additionally,  the  spectral  characteristics  of  reverberation  are 
essentially  the  sane  as  the  transmitted  signal,  the  intensity  of  reverberation 
varies  with  the  range  of  the  scatterers,  and,  finally,  the  intensity  varies 
with  the  intensity  of  the  transmitted  signal.  The  last  fact  ia  responsible 
for  the  lack  of  improvement  in  detectability  when  the  source  level  is 
increased;  i.e.,  relative  to  reverberation. 

The  literature  on  the  subject  of  scattering  and  reflection  of  sound  from 
boundaries  is  voluminous.  It  is  not  the  intention  of  this  report  to  survey 
the  entire  field.  There  are,  however,  some  key  references  which  form  the 
points  of  departure.  Fortuin1  has  surveyed  the  literature  up  to  the  beginning 
of  1970  on  reflection  and  scattering  of  sound  waves  at  the  sea  surface. 

Mackenzie,6  and  McKinney  and  Anderson7  have  considered  reflection  and 
backscattering  of  sound  from  the  ocean  bottom.  Their  work  is  a  cornerstone  of 
all  bottom  reverberation  work  to  date.  A  contemporary  collection  of  papers 
concerning  bottom  interacting  ocean  acoustics  appears  in  the  collection  edited 
by  Kuperman  and  Jenaen.6 

The  major  goals  of  this  survey  are  to  identify  the  methods  that  are 
currently  used  in  characterizing  reverberation  and  to  distill  the  essential 
mathematics  used  in  doing  so,  as  well  as  identifying  several  potential  new 
means  of  feiaracterizing  reverberation.  In  identifying  new  means  of  character¬ 
izing  reverberation,  attention  is  focused  on  obtaining  features  which  are  not 
used  in  the  mathematical  framework  in  the  past,  rather  than  identifying  only 
new  algorithms  or  approaches  to  make  the  standard  approaches  less  time  consum¬ 
ing  or  awre  comprehensive. 
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The  rewards  of  additional  aethods  or  features  by  which  one  can  charac¬ 
terize  reverberation  are  in  the  increased  ability  to  recognize  signals  or 
target  echoes,  as  well  as  obtaining  a  more  accurate  siaulation  of  reverbera¬ 
tion.  In  the  siaulation  arena,  there  are  two  aspects:  one  in  coaputer 
siaulation  for  testing  equipaent  performance  and  the  other  in  counteraeasure 
application;  e.g.,  generate  a  false  reverberation- like  field. 

Any  additional  methods  of  characterizing  reverberation  can  be  used  as 
predictors  of  other  reverberation  measures  as  well  as  in  casting  light  on  the 
aechanisms  that  generate  reverberation.  NOTE:  The  use  of  the  word  aeasure  in 
this  report  is  intended  for  the  mathematical  quantification  of  some 
attributes;  in  other  words,  a  metric  on  a  space,  distension  of  a  process  such 
as  topological  dimension,  a  statistic,  etc. 

The  prominent  descriptive  methods  in  reverberation  characterization  are 
given  in  a  group  of  papers  and  reports  that  originated  at  the  University  of 
Texas'  Applied  Research  Labs.2  “  10  11  14  11  14  “ 

A  review  of  reverberation,  scattering,  and  echo  structure  where  rever¬ 
beration  is  restricted  to  receiver  being  near  the  transmitter  is  given  by 
Claude  Horton,  Sr.9  Horton  deals  mostly  with  different  models  for  the 
scattering  from  rough  surfacea  and  points  out  the  assumptions  and  difficulties 
of  the  specific  models. 

The  reverberation  studies  used  for  characterizing  purposes  fall  into 
three  categories:  (1)  computer  simulation  from  which  reverberation  returns 
are  generated,  and  then  time  and  ensemble  averages  are  calculated;  (2)  model 
tank  observations  for  cast  surfaces  which  may  be  periodic  or  randomly  rough; 
and  (3)  observations  of  reverberation  at  sea  or  at  a  lake. 

An  example  of  computer  studies  is  the  simulation  of  reverberation  from  a 
moving  sea  surface  conducted  by  Bourianoff  and  Horton.16  The  sea  surface  is 
generated  by  using  a  100-term  Fourier  polynomial  with  random  phases.  The 
polynomial  is  of  the  form: 


100 

F(x,t)  =  I 
m=l 


2m+l 


C08(k2m+lX 


u>  t  + 
2m+l 


e2mH) 


E„  ,  are  the  random  phases  distributed  uniformly  in  the  interval  (0,  2n)  with 

ZDn’  1 

new  phases  generated  for  each  sample.  Amplitudes  are  generated  from  the  power 
spectrum  S(ui)  of  the  wave  spectrum  by: 
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for  • 


-Pierson1*  ware  spect 


«  _-p(w»  /•»)* 

S<»)  «jp  l*« 


t  *  980  cm/sec* 

*o  *  g/u  1/sec 

u  »  wind  velocity  la  cm/ sec 

a  *  8.1  z  l<f3 
o 

P  *  0.74  . 

The  results  of  that  simulation,  restricted  to  the  Neuman-Pierson  wave 
spectrum  and  pulses  which  are  short  compared  with  the  dominant  wave  length  of 
the  scattering  surface,  show  that  means,  variances,  and  covariances  are 
different  for  time  averages  as  compared  to  ensemble  averages.  In  the  study, 
it  was  seen  that  variances  computed  as  ensemble  averages  of  the  second  order 
moment  were  not  constant  and  hence,  the  reverberation  was  not  a  stationary 
process.  The  covariances  were  computed  for  the  tine  series  as  well  as  the 
ensembles  and  significant  differences  were  seen.  Specifically,  the  ensemble 
average  oscillated  more  rapidly  and  decayed  more  rapidly  than  the  time 
average  The  frequency  of  oscillation  of  the  individual  covariance  functions 
was  higher  than  the  center  frequency  of  the  transmitted  pulse.  The  behavior 
was  attributed  to  the  smallness  of  the  insonified  area  which  "produced"  a 
backscattering  coefficient  which  accentuated  the  higher  frequencies.  ’  This 
model  did  not  fit  experimental  results. 


STATISTICAL  THEORIES  OF  REVERBERATION  -  THE  MIDDLETON  APPROACH 

The  task  of  developing  a  deterministic  theory  of  reverberation  has  proven 
to  be  formidable.  There  is  the  theoretical  difficulties  of  solving  quite 
complex  boundary  value  problems  for  which  analytical  tools  are  only  partly 
developed.  There  is  also  the  practical  difficulties  of  identifying  all  the 
parameters  affecting  the  reverberation  process  and  obtaining  sufficiently 
accurate  measurements  of  these  parameters.  In  order  to  side  step  these 
problems,  Middleton  and  his  co-workers10  11  developed  a  semi-empirical  statis¬ 
tical  approach  which  will  be  discussed  below. 

In  characterizing  data,  two  papers  by  Plemons,  Shooter,  and  Middle- 
ton10  11  set  the  stage  for  analysis  of  signals  from  lake  and  sea  surfaces. 
These 'papers  will  be  referred  to  as  PSM. 

PSM  point  out  the  requirement  for  adequate  statistical  description  of  the 
reverberation  process.  There  are  two  observations,  really  assumptions ,  which 
have  some  experimental  support;  namely. 
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1.  The  reverberation  processes  are  generally  Gaussian  as  a  result 
of  large  number  of  scatterers. 

2.  The  covariance  function  provides  sufficient  description  of  the 
reverberation  processes  for  signal  detection  and  extraction  requirements. 

The  current  view  of  the  adequacy  of  these  assumptions  will  be  discussed 
later;  however,  these  form  the  starting  point  for  numerical  characterization. 
Another  key  element  discussed  by  PSM  is  the  need  for  suitable  ensemble 
averages.  Prior  investigations  typically  used  time-averages  on  single-time 
functions;  however,  since  the  dynamics  and  changing  environmental  conditions 
vary,  there  is  little  agreement  in  measurements  based  on  time  averages.  PSM 
claim  that  theirs  is  the  first  systematic  use  of  ensemble  averages  used  for 
detailed  analysis  of  non-stationary  mean  intensity  and  covariance  of  rever¬ 
beration.  The  second  order  statistics  are  used  to  study  and  test  Middleton's 
point  scattering  models.3 

The  point  scattering  models  of  Middleton  are  also  developed  by  Faure18 
and  Ol'shevskii. 19  One  may  collectively  call  it  the  FOM  approach.  The  FOM 
approach  is  based  on  the  essential  assumption  that  the  reverberation  is  the 
result  of  scattering  from  point  scatterers  located  independently  and  randomly 
in  space  and  re-radiating  incident  radiation  in  an  independent  manner  in  time. 
The  approach  is  quasi-phenomenological  in  that  there  is  no  a  priori  specifica¬ 
tion  of  the  impulse  function  of  the  point  scatterers;  hence,  no  a  priori 
calculation  of  reverberation  on  an  absolute  scale  can  be  performed.  Experi¬ 
mental  data  are  used  to  calibrate  such  models.  The  improvement  is  in  the 
elimination  of  the  need  for  boundary  conditions  and  solutions  of  reflection/ 
scattering  equations  with  stochastic  boundary  conditions. 

Since  the  statistics  that  Middleton's  models  yield  are  the  variances  and 
co-variances,  PSM  look  for  these  statistics.  Additionally,  the  data  sets  are 
validated.  The  validation  consists  of  tests  for  homogeneity  and  independence. 
The  validation  of  homogeneity  implies  the  stability  of  the  scattering 
mechanism;  in  other  words,  the  density  of  scatterers  is  examined;  fluctuations 
in  density  will  lead  to  data  that  are  inhomogeneous.  Tests  for  homogeneity 
include  the  Kolmogorov-Smirnov  two  sample  test  and  the  Wilcoxon  rank-sum  test. 
Other  tests  for  homogeneity  also  exist  and  will  be  re-examined  in  more  detail 
later. 

PSM  make  the  assumption  that  the  reradiative  properties  of  each  scatterer 
are  time  invariant  and  frequency  insensitive;  i.e.,  the  scatterers  are  assumed 
to  have  a  constant  dynamic  cross-section  which  depends  only  on  location 

and  direction  of  incoming  illumination.  The  assumption  also  includes  neglect 

of  Doppler  effects  due  to  platform  or  scatterer  motion.  The  scatterer  is 
represented  by  a  stochastic  filter  with  the  impulse  function 


V  v>  =  Y0lc^)  set  -  y 

where  y  .  is  the  dynamic  cross  section  for  scatterer  located  at  V(\,^).. 
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These  assumptions  lead  to  a  received  signal 

1 


U.(t)  *  S  h  (x/XJ  S(t  -  t)dt 
J  J  * 


% 


where  S(t)  is  the  signal  that  illuminates  scatterer  j  [due  to  the  6-function 
structure,  IK (t)  =  T0jS(t  -  A..)]. 

The  receiver  then  "sees"  the  reverberation  as  a  sum  of  all  the  basic  wave 
forms  expressed  as 


M 

X(t)  =  1  U(t.,x  ,0  ) 
j  J  J  J 


where  the  suaawtion  is  over  all  contributing  scatterers.  0.  is  the  set  of 
random  parameters  associated  with  j. 

How  if  reverberation  has  zero  mean  <X(t)>  =  0  for  all  t,  then  the 

6Q8 •  BV • 

complex  covariance  function  is  defined  by 


av. 


PSM  make  use  of  narrow  band  input  signals  which  are  represented  by  the 
complex  form 


S(\)  =  So(X)  exp(iwoX);  So(X)  =  Aq{\)  e^(M 

So(X)  is  the  cosqilex  enveloped  of  S(\) 

Aq(\)  is  the  real  envelope 
♦(A)  is  the  real  phase. 


A^(X)  and  #(X)  vary  slowly  with  respect  to  eiuJo\ 

The  quantity  Hq(t)  is  the  complex  envelope  of  the  covariance  function  and 
is  defined  as 


S  ” 


•  1 
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Ho(T)  5  S  So(X)*  So(X  +  x>  dX  5  «0W  expti^Ct)] 

•OO 


K0  =  1/2 |Ho(l) I 

♦c  =  -t*Q'l«»H0/lteH0) 

thus 


are  real  quantities 


K(x/tj)  =  GCtj)  Kq(t)  cos[woX  +  ♦(X)] 


where  G  is  a  function  which  is  not  determined  in  PSM  but  where  t  ormalized 
covariance  function  extracts  and  eliminates  the  need  for  this  funct 

The  normalized  covariance  function  is  then 


<K>n(x  |t1)  =  [2<K>o(x)/Ej  cos(woX  +  <*>Q(x)) 


00 

E  is  the  energy  in  the  complex  transmitted  signal  =  /  |S(A)  |2d\. 

-00 


Another  requirement  is  the  prevention  of  intermixing  reverberation 
returns  which  is  achieved  by  proper  spacing  of  output  signals  or  pings.  This 
requirement  carries  the  label  of  independence  of  observations. 

The  demonstration  of  homogeneity,  i.e.,  that  data  come  from  the  sasie 
parent  population,  and  independence,  i.e.,  that  data  contain  no  contamination 
of  one  observation  by  another,  is  required. 

PSM  and  other  ARL  investigators  use  the  quadrature  sampling  techniques  to 
convert  analog  signals  to  digital  signals.  This  technique  was  developed  by 
Grace  &  Pitt.12  The  reverberation  return  is  expected  to  be  narrow  band  since 
the  transmitted  signal  is  narrow  band  and  thus  the  reverberation  return  is 
expressed  as: 


X(t)  =  X  (t)  cos  u>  t  +  X  (t)  sin  w  t 
c  os  o 

u>o  =  transmit  center  frequency 

X  ,X  =  quadrature  components  which  vary  slowly  compared  with 
cos  U)  t. 
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Ensembles  at  time  t^  are  constructed  by  selecting  the  value  of  some 

attribute,  e.g.,  amplitude,  at  a  time  t^  from  each  reverberation  event. 

Different  ensembles  are  formed  by  choosing  different  ti* 

An  ensemble  statistic  may  then  be  constructed;  for  instance,  the  estimate 
of  a  mean  amplitude  at  time  t^  would  be  evaluated 

1  N 

<X(t,)>  ens.  av.  =  ^  I  X.(t.) 

1  "  i=l  1  1 

where  N  is  the  number  of  reverberation  observations  or  records  that  are  taken. 
An  experimental  covariance  function  can  be  evaluated  as 


<K>  (trt2)  =  J  I  X.Ctj)  X.(t2). 
i=l 

If  one  follows  PSM's  approach  for  narrow  band  signals,  the  covariance 
function  may  be  expressed  as 

<K>  (tj,t2)  =  ClO^tj.tj  +  X)  cosier  +  «|»o  (tj.tj  +  X)] 

where 


<V  (ti»V  =  I<R>x  (ti»V  +  <A>x 


where  <R>X  and  <A> ,  the  quadrature  components  of  the  covariance  are: 


<R>  = 


1  1 


N 


«  ■  5  •  S  £  lXc,l<tl>  xc,i(t2>  + 

«  *  5  •  S  jj  '-Xc.i(t2)  *  xc,i(t1)  Xs,i(t2» 


leading  to  the  estimate  <$>  of  the  phase  0  according  to 
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Local  statioaarity  is  tested  via  the  requirement  that 


+  x)  =  ■  x) 


hold;  i.e.,  compare  the  phase  for  positive  values  of  x  with  the  negative  image 
of  the  phase  when  X  is  negative. 


EFFECT  OF  ENSEMBLE  SIZE  ON  COVARIANCE 

PSM  pay  attention  to  the  fact  that  experimentally  one  deals  with  a  finite 
number  of  observations  and  ensemble  averaging  is  thus  performed  over  a  finite 
number  of  reverberation  returns .  The  experimental  co-variance  is  therefore 
only  an  estimate  of  the  true  infinite  sized  ensemble  co-variance.  The  vari¬ 
ance  of  the  co-variance  should  decrease  as  0  (M  *)  where  M  is  the  size  of  the 
sample  ensemble.  PSM  argue  that  a  sample  size  of  150  members  is  acceptable. 

PSM  also  study  the  estimate  of  covariance  using  time  averages  where  the 

time  average  over  a  single  member  must  be  used.  For  the  ensemble  member , 
the  time  averaged  co-variance  then  appears  as 


CjCxlTj,^)  =  |  f  X  it)  X.it  +  X)  dt;  T  =  T2  -  11 
T1 


in  order  for  c^CxlT^T^)  to  equal  the  ensemble  covariance  the  following 
requirements  need  to  be  met:  (1)  T  -*•  •,  (2)  stationarity,  and  (3)  ergodicity. 

Other  works  to  originate  from  ARL  concerning  the  subject  of  reverberation 
and  in  the  line  of  the  Plemons,  Shooter,  and  Middleton  approach  include  those 
of  Frazer13  and  Wilson.2  The  objective  of  Frazer's  study  is  the  attainment  of 
a  more  precise  underlying  statistical  distribution  of  the  scattered  fields. 
Frazer  points  out  that  the  majority  of  studies  focused  on  the  examination  of 
surface  backscattering  strength  with  assumptions  that,  as  a  result  of  large 
number  of  scatterers,  the  outputs  of  individual  receiver  elements  should  be 
zero  mean  Gaussian  distributed.  The  conclusions  of  his  studies  include  the 
desire  for  a  stronger  delineation  of  the  conditions  for  producing  non-Gaussian 
reverberation  and  detecting  it. 

Wilson's  study2  continues  along  the  lines  of  ARL's  earlier  work  in 
verifying  and  calibrating  the  Middleton  model  of  reverberation.  Expanded 
experimental  attention  is  given  to  the  study  of  horizontal  and  vertical  co- 
variance;  the  first  four  univariate  moments  of  the  reverberation  are  examined 
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and  normality  tests  are  applied.  Additionally,  measurement  and  interpretation 
of  multiple  receiver  reverberation  are  investigated.  Specific  results  on  the 
covariance  include  (1)  decrease  of  a  value  of  0.1  within  a  separation  of  four 
wave  lengths  of  the  normalized  envelope  o'  the  horizontal  covariance;  (2)  the 
lower  bound  of  the  vertical  covariance  as  •?  function  of  spectral  separation  is 
0.2;  (3)  the  first  four  moments  of  mean,  variance,  skewness,  and  kurtosis 
contained  an  oscillatory  component  which  depended  on  the  order  of  the  moment; 

i.e.,  the  n^  moment  appeared  to  oscillate  at  fundamental  frequencies  of  nu>o 
where  u>q  is  the  transmit  frequency;  and  (4)  ensembles  were  found  to  be  non¬ 
normal  due  to  slight  skew  and  large  kurtosis.  Other  details  of  the  structure 
of  the  covariance  envelopes  and  phases  are  presented  which  reinforce  the  claim 
of  dependence  of  amplitude  and  co-variance  structure  on  the  geometries  of 
transmitter/ receivers,  pulse  widths,  frequencies,  and  other  parameters. 

It  is  appropriate  to  note  here  another  fundamental  investigation  which 
appeared  in  report  form  only;  namely,  the  work  of  J.  Blue.14  Blue's  report  in 
effect  set  the  stage  for  the  other  ARL  works  of  the  1970s.  It  is  more  opera¬ 
tionally  oriented  since  it  culminated  in  the  evaluation  of  the  performance  of 
signal  processors  against  the  recorded  reverberation. 

Blue's  report  is  of  specific  importance  since  it  obtains  reverberation 
generated  by  bottom  backscattering.  Most  of  the  reverberation  work  is  for 
surface  generated  backscattering.  The  signals  used  by  Blue  are  1  msec, 
80  kHz,  cw  pulses,  where  the  projector  and  receiver  arrays  are  positioned 
9  feet  off  the  bottom,  for  a  3-degree  beam  width  at  -3  dB.  Blue  uses  the 
Edgeworth  series  to  describe  departures  from  Gaussianity  of  the  amplitude 
distribution,  or  departures  of  the  envelope  distribution,  from  the  Rayleigh 
distribution.  Reasons  for  such  departures  are  the  insufficiency  of  ensemble 
sizes  and  the  possibility  of  a  small  number  of  scatterers  being  involved  in 
scattering  at  any  one  time.  These  conjectures  are  supported  by  the  observa¬ 
tions  that  the  amplitudes  departed  from  the  Gaussian  distribution  and  corres¬ 
pondingly  the  envelopes  were  not  Rayleigh  distributed.  Blue  attributed  the 
non-Gaussian  features  of  the  "operational”  ensemble  data  to  the  inhomogeneous 
character  of  the  returns.  A  normalization  procedure  was  applied  to  the  raw  or 
"operational"  data.  The  instantaneous  values  were  normalized  over  3  msec 
locally  stationary  portions  of  each  return  by  dividing  by  the  RMS  estimates  of 
each  return.  This  procedure  led  to  good  fits  with  the  Gaussian  distribution. 

The  comparisons  of  square  law  and  multiplicative  detectors  such  as 
inhomogeneous  reverberation  demonstrated  that  while  the  square  law  detector 
was  superior  to  the  multiplicative  detector,  the  detection  loss  is  lower  for 
the  multiplicative  detector.  Blue  also  shows  that  inhomogeneous  reverberation 
degrades  the  performance  levels  obtained  in  homogeneous  reverberation. 

Another  investigation  of  bottom  reverberation  is  that  of  Pitt.ls  Pitt's 
results  show  that  a  smooth  estimate  for  the  backscattering  was  difficult  to 
obtain,  even  with  good  physical  conditions.  Pitt  stresses  that  the  ensemble 
statistics  obtained  were  independent  of  the  area  sensed  for  a  wide  range  of 
areas.  The  short  cw  and  long  linear  FM  pulses  appear  to  be  distributed  with 
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the  same  Rayleigh-like  forms.  The  descriptive  tools  were  plots  of  relative 
level  (-dB )  versus  grazing  angles,  estimated  backscatter  coefficient  versus 
grazing  angles,  envelope  statistics  of  number  of  occurrences  versus  asq>litude 
in  form  of  histograms  fitted  with  Rayleigh  distribution,  and  output  S/R  levels 
versus  input  signal/reverberation  (S/R). 

The  paper  of  Lord  and  Plemons20  undertakes  the  task  of  characterizing  and 

simulating  signals  reflected  from  the  sea  surface  in  the  forward  direction. 

The  sea  surface  is  characterized  by  a  roughness  parameter  g  =  4n6  sin  B/X; 

6  =  rms  surface  wave  height,  X  the  acoustic  wave  length,  and  6  the  grazing 

angle.  The  signals  are  quadrature  sampled.  A  random  variable,  £  =  E/E  is 

used  to  characterize  returns  where  E  =  energy  of  surface  scattered  pulse,  and 

E  is  the  energy  corresponding  to  scattering  or  specular  reflection  from  a 
sp 

flat  surface.  E  serves  as  an  upper  bound  or  normalization  constant, 
sp 

A  unimodal,  hypothetical,  probability  function  is  uqed  to  simulate  the 
distribution  according  to 


f(4;a,P)  =  iW(P)  t*'1  (1  "  °P1  for  0  -  1 1  1 


where  f  =  0  elsewhere.  The  parameters  in  this  theory  are  then  a  and  p.  The 
means  and  variance  of  £  can  be  calculated  for  the  Beta  distribution  (the 
function  given  before)  and  are: 


M£  =  «/(<*  +  0) 

6^  =  op/(o  +  P)2(or  +  p  +  1)  . 


The  parameter  estimates  a  and  p  are  then  obtainable  as: 

a  =  P£  Km^Ict^MI  -  M^)  -  1] 


P  =  (1  -  t  CM^ I ) (1  ”  Mj)  -  H 


where  \xr  and  af 


are  the  sample  estimates. 


Lord  and  Plemons  then  compare  the  observed  data  and  postulated  Beta  dis¬ 
tribution  via  the  Kolmogrov-Smirnov  test  for  the  cumulative  distribution 
functions 
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A  A  ^  -A  A 

cdf  of  f  *  F(£;a,0)  *  J  f(x;a,p)  dx 


The  empirical  cdf  is  computed  from  the  ordered  sasples.  The  K-S  test  indi¬ 
cated  that  there  was  no  reason  to  reject  the  hypothesis  of  a  Beta  distribution 
at  the  5  or  10  percent  levels  of  significance.  The  prise  use  of  this  result 
is  claimed  to  be  in  Monte  Carlo  siaulation. 

Additional  structure  in  the  pulses  is  investigated  via  the  assumption  of 
autoregressive  or  recursive  features.  The  assumption  is  supported  by  the 
apparent  temporal  coherence  of  the  forward  scattered  pulse. 

The  linear  autoregressive  model  used  is 


Sr(j,D  =  Z(j,l) 

*k 

SCj.k)  *  2  A(k,£)  S  (j ,k  -  A)  +  Z(j,k) 
r  £=1 

E£>(  ,  )  represents  surface  reflected  signals 

Z(  ,  )  represents  zero-swan  additive  random  variables 

A(  ,  )  is  the  set  of  regression  constants. 

j  ranges  over  the  number  of  pulses  in  the  ensemble  1  *npt  n^  is  the  time 

increment,  and  n  is  the  number  of  pulses  (1  < 

P 

order  of  the  autoregressive  process. 

The  least  squares  estimators  are  obtained 


j  <  n  ,  2  £  k  £  n  ).  1  is  the 

P  t  K 

by  minimizing  the  squared  error 


S2(k)  =  2 
j=l 


lSr(j,k)  - 


2 

A*  1 


A(k,£)  Sr(j,k  -  A) I* 


where  -successive  sets  of  equations  are  formed  by  setting 


9S2(k) 

8A(k,n) 


0  . 
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Lord  and  Plemons  demonstrate  with  this  scheme  that  the  quadrature  compo¬ 
nents  at  any  given  time  are  largely  determined  by  their  immediate  predecessors 
and  less  so  on  the  distance  past.  Thus,  while  recognizing  that  no  insight 
into  the  scattering  mechanism  is  obtained,  it  is  shown  that  a  first  or  second 
order  linear  autoregressive  model  provides  a  reasonable  simulation  of  the 
non-stationary  surface  reflected  signals. 

Other  examples  of  reflection  and  scattering  models  and  experimental 
results  are  available;  however,  it  is  not  the  purpose  here  to  review  material 
that  is  already  reviewed  by  Fortuin1  and  Horton.9  It  is,  however,  appropriate 
to  note  the  existence  of  the  Marsh-Kuo  approach  given  in  some  papers  in  the 
early  60s.21'22  23 > 24 


GEOACOUSTIC  MODELS  AND  ECHO  CORRELATION  TECHNIQUE 

The  statistical  approch  of  Middleton  described  in  the  last  section  is 
complemented  by  the  geoacoustic  models,  and  echo-echo  correlation  techniques 
possible  with  aid  of  contemporary  technology.  The  geoacoustic  model  used  by 
Russian  investigators  described  here  is  useful  in  understanding  the  method¬ 
ology  required  for  incorporation  of  multiple  scales  of  inhomogeneities. 

Foreign  investigations  of  sound  scattering  by  the  ocean  bottom  follows 
the  experimental  investigations  cited  before.  The  review  of  Bunchuk  and 
Zhitkovskii25  is  representative  of  available  foreign  results.  The  foreign 
papers  do  not  offer  any  novel  information  or  methodology;  however,  it  is 
illuminating  to  observe  the  similarities  in  approaches. 

Ivakin  and  Lysanov28  investigate  bottom  reverberation  and  develop27  a 
geoacoustic  model  of  the  bottom  with  a  view  of  verifying  the  insensitivity  of 
the  results  of  Reference  27  to  varying  the  form  of  the  correlation  coefficient 
of  the  refractive  index  fluctuations. 

The  geoacoustic  model  assumes  a  bottom  thickness,  h.  The  sediment  is 
regarded  as  a  water  saturated  medium  where  the  average  values  of  sound 
velocity  and  density  differ  only  slightly  from  one  sublayer  to  another  and  the 
upper  sedimentary  layer  differs  very  slightly  from  the  water-only  layer.  The 
result  of  this  model  is  that  there  is  no  reflection  from  the  upper  sedimentary 
interface  at  all,  but  rather  scattering  appears  from  inhomogeneities  in  the 
sediment  layer. 

The  coefficient  of  backscattering  from  an  absorbing  layer  with  random 
inhomogeneities  is 


Mg  =  (2P)"1  Mysin  «|f 


where 


P  =  sound  absorption  coefficient 
M^  =  volume  scattering  coefficient  in  the  non-absorbing  medium. 
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The  above  expression  is  valid  for  20h/sin  X  »1.  If  the  linear  dimensions  of 
the  scattering  volume  are  larger  than  the  inhosiogeneity  scale,  the  volume 
scattering  can  be  related  to  energy  spectrum  G(q)  by 


H  =  2nk4  G(q) 


where  the  energy  spectrum  is  given  by 


G(q)  =  <M2>  /  N(r)eiq*r  dr 


where  <p*>  and  N(r)  are  the  mean-square  value  and  correlation  coefficients  of 
the  refractive  index  fluctuations;  q  =  £  -  kQ,  k,kQ  are  the  incident  and 
scattered  wave  vectors. 

Observations  indicate  that  the  abaorption  coefficient  is  proportional  to 
frequency;  thus,  for  frequency  independent  scattering  coefficient  the  energy 

.3 

spectrum  must  be  proportional  to  k  Anisotropic  inhomogeneities  are  of  a 
large  acale  in  the  horizontal  direction  and  small  in  the  vertical-depth  scale. 
The  correlation  coefficient  is  then  decomposed  as 


which  leads  to 


N(r)  *  MjCz)  •  Ha(p)  p  =  (x,y) 


G(q)  *  <M*>  Gj(y)  G2(v) 


where  G^(y)  is  the  normalized  energy  spectmm  of  small  scale  inhomogeneities, 
G2(v)  of  the  large  scale  inhomogeneities. 

For  example: 


°i(»>  ■  h  i  *i(,)  *lYi  *  s  2  "i<,> 


0*<V)  *  J  F.(p)  eiVp  dp  . 
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G1(y)  is  essentially  a  constant  due  to  small  scale  assusiption  leading  to 
eiYz  S  1. 

For  horizontally  isotropic  inhomogeneities,  the  following  form  is  used: 


N2(p)  =  exp(-  |p  |/pQ)  p  =  V*z  +  yz 

which  yields 


g2(v)  =  Jr  a  +  v2Po>'3/2 


and  for  vpQ  »1,  G2(v)  =  (2npov3)”  . 

The  backscattering  case  v  =  2k  cos«J»  leads  to 


G2(v)  =  l/(167ipok®  cos3t|0 


Mg  =  <p2>(k^)(z/PQ)  ®i“  «J</l6ncoss»|i. 


One  is  now  interested  in  varying  the  form  of  the  correlation  coefficient. 
Operationally,  this  is  justified  since  different  bottoms  will  be  described  by 
different  conditions.  Ivakin  and  Lysanov  recommend  the  use  of  a  certain  class 
of  functions  wherein  a  change  of  form  of  the  correlation  coefficient  does  not 
alter  the  fundamental  results.  The  class  of  function  investigated  is  that 
where  N2(p)  has  a  non-vanishing  first  derivative  for  p  =  0.  This  class  can  be 

~  -3 

shown  to  maintain  an  energy  spectrum  G2 (v) ~v  ,vpo»l.  Since  in  reality  there 

are  no  infinite  gradients,  in  addition  to  the  inhomogeneity  scale  pQ,  an 

additional  parasieter  px  characterizes  the  minimum  thickness  of  the  transition 
layers  between  individual  inhomogeneities  and  has  dimensions  of  length. 

V<M>a/Pi  determines  the  limiting  gradient  of  the  random  field  p  in  the  medium. 
The  function 


N2(p)  =  (e‘lpl/po  -  ee"lp|/pl)(l  -  e)'1,  e  h  Pl/pr 


will  have  a  zero  derivative  for  p  =  0  and  differ  from  single-scale  correlation 
only  in  the  vicinity  of  zero. 


5 


NCSC  TM  362-82 


The  energy  spectrum  that  corresponds  to  the  above  function  is 

c  -  P°  r _ 1 _  “  _ e2  , 

2  2n(l  -  e)  (1  +  v2pQ2)3/2  (1  +  v2p2)3/2 


In  the  interval  l/pQ  <  V  <  l/p  ,  then  it  can  be  shown  that  the  models  yield 

invariant  characterization  of  the  backscattering.  These  results  are  further 
extended  to  multiple  scale  models. 

A  somewhat  different  approach  altogether  is  taken  by  Dunsiger,  Cochrane, 
and  Vetter.2*  Dunsiger  et  al  were  interested  in  characterizing  the  seabed 
using  ping-to-ping  echo  coherence.  The  model  they  used  is  a  two  process  model 
where  a  deterministic  but  unknown  reflection  filter  operates  in  parallel  with 
a  random  scattering  filter.  They  remove  large  scale  topography  effects  by 
realigning  the  echoes  to  a  standard  time  reference.  Then  coherence  versus 
frequency  is  studied  using  a  broad  spectrum  impulsive  excitation  which,  after 
processing,  yields  data  in  bandwidth  from  less  than  1  kHz  to  more  than  10  kHz. 
The  coherence  or  correlation  approach  yields  reasonable  results  for  soft 
bottoms  while  harder  bottoms  yield  relations  which  are  inconsistent  with  any 
master  curves  derivable  from  a  linear,  parabolic,  or  Gaussian  correlation 
function. 

Other  studies  of  the  dependence  of  spatial  and  temporal  correlations  of 
scattered  underwater  sound  are  those  of  Clay  and  Medwin29  30  who  study  the 
correlation  for  forward  scattering  when  source  and  receiver  are  far  away  and 
surface  slopes  are  assuswd  to  be  small  such  that  no  shadowing  effects  occur. 
The  experimental  part  was  conducted  in  a  test  tank  at  the  Haval  Postgraduate 
School . 

It  is  thus  quite  apparent  that  a  substantial  amount  of  theoretical  devel- 
opment  and  experiments  designed  to  verify  them  are  based  on  correlation  inves¬ 
tigations  where  coherences  or  covariances  are  measured  either  for  amplitudes 
alone  or  amplitudes  and  phases.  General  studies  on  covariances  coherences  and 
spectral  analysis  from  the  pure  mathesuitical  point  of  view  are  given  by 
Middleton31  and  Nuttal.32  Middleton's  book31  contains  a  thorough  discussion 
of  many  aspects  of  statistics,  probability,  and  coherence/covariance  analyses 
as  they  apply  to  communications  systems.  These  methods  form  the  basis  of  the 
work  of  ARL  described  earlier. 

In  this  context,  it  is  appropriate  to  note  the  report  of  Green33  who 
investigated  the  target  signature  of  a  mine  for  large  bandwidth-time  (BT) 
product  with  greater  than  CW  pulses.  The  replies  correlator  is  accepted  as 
the  optimum  filter  for  an  echo  embedded  in  Gaussian  white  noise  and  when  the 
echo  exactly  replicates  the  transmitted  pulse.  The  Echo-Echo  Correlation 
technique,  as  the  name  implies,  attempts  to  correlate  one  echo  with  another 
and  is  particularly  adapted  to  account  for  experimental  and  target  induced 
variabilities . 
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Green  emphasizes  that  there  is  a  great  deal  of  variability  in  both  the 
echo  shape  and  echo  power.  The  attempt  to  obtain  more  target  structure  by 
using  higher  resolution  pulses  did  not  succeed  as  projected.  It  was  observed 
that  for  a  given  pulse  bandwidth,  Echo-Echo  Correlation  requires  three  times 
the  same  BT  product;  also,  EEC  does  not  provide  precise  range  information 
available  with  RC.  However,  EEC  is  claimed  to  be  advantageous  over  RC  for  the 
following  capabilities:  (1)  simple  processing  using  time  domain  correlation 
in  real  time,  (2)  reduced  display  requirements,  (3)  stable  correlation  perfor¬ 
mance  with  increased  bandwidth  (RC  does  not),  (4)  improved  false  alarm  rate, 
(5)  better  motion  tolerance,  and  (6)  better  suited  for  use  with  Pseudo  Random 
Noise  (PRN)  signals. 


LEVEL  CROSSINGS  ANALYSIS 

A  completely  different  approach  to  study  of  reverberation  is  possible 
with  the  use  of  zero-crossing  and  level-crossing  analysis  techniques.  The 
most  important  feature  of  level-crossing  analysis  is  that  it  is  iaiplicitly 
focused  on  features  relating  to  phase  rather  than  amplitude.  The  zero 
crossing  approach  is  that  where  the  rate  of  occurrence  of  zero-level  amplitude 
for  noise  and  signal  •  noise  are  compared.  Rice34  investigated  the  distribu¬ 
tion  of  crossings  of  the  y  =  0  axis  of  a  stationary  Gaussian  noise  function 
y(t)  and  derived  the  following  formula  for  the  average  number  of  zero  cross¬ 
ings  per  unit  time 


A  =  21  J  f%(f)  df/  J  uj(f )  df]5* 
o  o 


where  u>(f)  is  the  spectral  density  of  the  function  y(t).  Bom  and  Conoly3S 
extended  this  work  to  develop  the  zero-crossing  shift  as  a  detection  method. 
Their  work  had  been  commented  upon  by  Riter  and  Boatright36  and  replied  to  by 
Bom  and  Conoly.37  They  recognize  the  importance  of  investigating  the  distri¬ 
bution  of  time  intervals  x  between  consecutive  zero  level  crossings  and  point 
out  Rice's  observation  of  the  substantial  difficulty  involved.  Bom  and  Conoly 
claim  that  at  the  time  of  their  work  no  one  had  found  a  satisfactory  theoreti¬ 
cal  solution,  and  hence  they  made  use  of  an  empirical  approach  for  finding  the 
probability  density  function  p(x)  that  would  approximate  the  experimentally 
observed  distribution. 

The  procedure  used  by  Bom  and  Conoly  determined  the  zero  crossing  by 
passing  a  noise  signal  and  counting  the  number  of  occurrences  within  a 
selected  range  t  +  At  over  a  fixed  time.  The  probability  density  of  Pearson 
type  III  family  was  used  for  convenience: 


p(t)  =  [p/r(m)]  (px)"_1e'pX 


with  p  >  0  and  m  >  0  constant. 
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This  type  of  distribution  coupled  with  an  assumption  of  very  rapid  decor¬ 
relation  between  successive  intervals  makes  the  zero  crossing  a  renewal 
process . 

For  ■  =  1 ,  p(x)  =  pe  and  the  probability  that  n  more  crossings  occur 
before  t  =  T  is 


Ufl(T)  =  [ (pT)n/n! ]e_pT 


as  pT  ■*  the  estimates  of  crossing  rate  n/T  are  asymptotically  normally  dis¬ 
tributed  about  p  with  variance  p/T.  In  the  general  case  p  =  p(p,m)  above,  we 
have  an  asymptotically  normal  distribution  about  p/m  with  variance  p/m2T.  A 
series  of  1-second  counts  is  then  made  £*,  £2 ,  £3...£^  and  the  mean  rate  eval¬ 
uated  as 


For  a  large  enough  T,  the  estimate  Z  forms  a  normal  distribution  about  a  mean 
value  A  =  p/m  and  then  Z  can  be  considered  as  independently  drawn  from  a  dis¬ 
tribution  N(A,a).  When  the  spectral  content  changes  slowly  so  that  the  dis¬ 
tribution  can  be  assumed  stationary  over  the  sampling  interval,  it  is  found 
that  the  mean  changes  faster  than  the  variance  of  the  distribution;  and  the 
resulting  distribution  can  be  regarded  as  the  initial  distribution  N(Ao>o) 

with  a  shifted  seas;  i.e.,  N(Alto)  and  a  Neyman-Pierson  detection  scheme  can 
be  used. 

Riter  and  Boatright  comments  generalize  to  discussion  of  level-crossing 
detectors.  The  detectors  count  the  rate  at  which  a  signal  crosses  a  given 
threshold  in  the  positive  direction  or  up-going  direction.  They  refer  to 
Papoulis3*  who  shows  that  the  level  crossing  rate  A^Ct)  is 


Aa(T)  *  *  VT) 


e-a*/2R(o) 


where  A  (x)  is  the  number  of  zero  crossings  in  time  x  and  R(o)  the  autocorre- 
o 

lation  of  the  signal  •  noise  at  t  =  0. 

The  concept  of  level  crossing  also  appears  in  an  informal  report  authored 
by  McLeroy,  Leadbetter,  and  Wegman.39  The  study  uses  two  statistical  tech¬ 
niques;  namely,  comparisons  of  zeores  of  reflected  waveforms  and  spectral 
comparisons  of  the  transmitted  and  reflected  waveforms.  The  targets  used  were 


NCSC  TM  362-82 


a  small  sphere  and  a  larger  sheet  metal  object  formed  of  short  concentric 
cylinders  stationed  on  a  single  axis.  The  zero  pattern  crossings  were 
summarized  by  means  of  histograms  of  length  intervals  between  successive 
zeroes.  The  zero-crossing  methods  involved  the  concepts  of  instantaneous  wave 
and  frequency.  An  arbitrary  waveform  x(t)  is  represented  by  the  real  part 

A  A 

in-phase  component  of  the  waveform  x(t)  +  ix(t)  where  x(t)  is  the  quadrature 
component  of  x(t)  (Hilbert  transform).  One  can  use  the  representation  in 
terms  of  a  principal  value  integral 


ds 


The  instantaneous  phase  i|>  is  then  the  argument  of  [x(t)  +  i  x(t)]  and  the 
instantaneous  frequency  is  d/dt  ijj(t). 

A 

Typical  examples  are  if  x(t)  =  cos  tut,  then  x  (t)  =  sin  u>t  (the  quadra¬ 
ture  component. 

Instantaneous  phase  =  avg  (cos  tut  +  isin  tut)  =  avg  (e1U,t)  =  tut. 
Instantaneous  frequency  =  d/dt  (utt)  =  tu. 

One  may  represent 


x(t)  =  J  x+(A)  e^^dA. 


~  o  . .  oo 

x(t)  =  i  J"  x+(A)  e1  tdA-i  J  x+(A)  e1  td\ 


x(t)  +  i  x(t)  =2  J  x+(\)  e^^tdA 


tj>(t)  =  Avg  {2  J  x+(\)  ei^td\} 
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We  nay  then  write  x(t)  +  i  x(t)  =  A(t)  e*^*^ . 
could  observe  that  ty(t)  =  2nN(t)  with  N(t)  being  the 

a 

the  origin  by  the  complex  waveform  x  +  i  x  in  (o,t). 
not  backtrack,  then  N(t)  would  be  the  number  of  zero 
phase-space  diagram  here. 

A 

X 


X 


Other  references  to  the  zero  crossing  concepts  include  the  phase  filter  zero 
crossing  counter  of  Huggins  and  Middleton40  41  and  analysis  of  a  narrow  band¬ 
pass  filter  of  Bendat.42 

Bendat42  calculates  the  average  number  of  zero  crossings  per  second,  at 
the  output  of  a  narrow  bandpass  filter  of  rectangular  shape,  when  there  is  a 
sine  wave  in  Gaussian  noise,  to  be 


If  we  take  x(o)  =  o,  one 
number  of  revolutions  of 

A 

If  the  vector  (x,x)  does 
crossings.  Refer  to  the 


nQ  =  2fo{[S/N  +  1  +  (ffi2/12fo2)]/[S/N  +  \)}h 


where  S/N  is  the  signal-to-noise  power  ratio 
fg  the  filter  bandwidth 
fQ  center  frequency  of  the  filter. 

There  are  other  zero-crossing  information  displays.  The  periodmeter43 
measures  the  period  of  each  cycle,  i.e.,  the  distance  between  every  other  zero 
crossing,  and  is  displayed  as  a  height  off  an  arbitrary  base  line.  When  noise 
alone  is  present,  a  random  distance  distribution  from  base  line  is  present 
while  for  a  signal  the  distance  between  zero  crossing  is  more  uniform  and  the 
scatter  of  points  is  reduced. 

In  discussing  the  nature  and  characterization  of  reverberation,  it 
becomes  apparent  that  the  characterization  is  dependent  upon  the  detection 
concepts  or  at  least  that  the  detector  used  for  detecting  signals  is  crucial 
in  determining  what  features  of  reverberation  should  be  studied.  One  impor¬ 
tant  concept  in  detection  is  that  of  Robust  Detection,  which  is  used  by 
Dwyer,44  Poor  et  al,45  and  Martin  and  Schwartz.46 

The  robust  detector  performs  well  over  a  range  of  noise  or  reverberation 
distributions.  Typically,  a  robust  detector  will  outperform  a  detector  which 
uses  an  optimum  Gaussian  procedure  when  the  underlying  distribution  of  noise 
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or  reverberation  is  non-normal.  When  the  underlying  distributions  are  normal, 
a  robust  detection  procedure  performs  almost  as  well  as  an  optimum  Gaussian 
procedure.  Robustness  is  typically  evaluated  by  considering  the  asymptotic 
relative  efficiency  (ARE)  versus  £  where  e  is  a  perturbation  parameter  in  the 
cumulative  distribution  function  of  the  noise  of  reverberation  process.  The 
salient  feature  of  the  ARE  is  that  it  is  a  measure  of  the  local  asymptotic 
performance. 46 

Earlier  in  this  review,  points  that  merit  detailed  discussion  were  iden¬ 
tified.  These  are  (1)  adequacy  of  the  assumption  of  Plemons,  Shooter,  and 
Middlteon3  11  concerning  the  Gaussianity  assumption  and  the  sufficiency  of  the 
covariance  function  and  (2)  other  statistical  tests  for  features  such  as 
homogeneity  of  data.  These  topics  will  be  explored  in  greater  detail  in 
future  reports  dealing  with  new  models  and  new  reverberation  characterization 
techniques.  Special  attention  will  also  be  given  to  the  topic  of  zero  cross¬ 
ings  and  level  crossings  making  substantial  use  of  the  work  of  Cramer  and 
Leadbetter47  and  subsequent  work  by  Leadbetter. 
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